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Abstract 

This work presents fluid dynamics videos obtained via numeri- 
cal (CFD) calculations using ComSol (finite elements method) soft- 
ware, showing the evaporation of HFE7100 (3M company refriger- 
ant) into a nitrogen gas flow along the liquid interface. The overall 
temperature evolution and liquid motion, which is caused by surface- 
tension (Marangoni) and buoyancy (Rayleigh) instability mechanisms, 
are shown as well. Flow behavior in the liquid caused by the afore- 
mentioned instability mechanisms can be nicely seen. Finally, these 
observations are made for three liquid thicknesses in order to appreci- 
ate the qualitative influence of confinement. 

1 Introduction 

This work falls into a general framework which consists of observing the 
behavior of patterns and structures that can be formed after instability 
onset in an evaporating liquid layer. In previous work, we studied theoretical 
instability thresholds in pure fluids [1,2] and in binary mixtures [3,4]. What 
is of interest here, is a two-dimensional numerical simulation study of the 
transient temperature and fluid motion in the liquid for a liquid evaporating 
into a nitrogen gas flow. The chosen liquid is HFE7100 (an electronic liquid 
produced by 3M). The numerical (CFD) simulations are performed using 
the software ComSol (finite elements method). The evaporation causes the 
instability and the gas flow evacuates the liquid vapor. The setup used for 
this numerical simulation is represented in Fig. [T] and is inspired from the 
CIMEX experimental setup of ESA [5]. 
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Figure 1: A scheme of the system 



The gas flow is maintained at 100 ml/min in a channel of 3 mm height, 
while three different liquid thicknesses are considered: 2, 4 and 8 mm. The 
width of the whole setup is 50 mm. The cover between the liquid and gas 
channel is 200/xra thick. At the middle of this cover, there is an opening with 
a width of 10.6 mm, allowing contact between the liquid and gas channel. 
These items are taken into the geometry of the numerical software ComSol. 
The boundaries of the whole system are kept at an ambient temperature 
and pressure of respectively 298 K and 1 atm, except for the gas channel 
outlet where only the ambient pressure is respected. Also, the whole system 
is surrounded by walls except for the gas flow inlet and outlet. The inter- 
face is kept at a constant height, since in the ESA experimental setup the 
liquid is to be replenished at the same rate as the evaporation rate. At the 
interface, flux conservation is maintained and a tangential stress balance is 
considered. Furthermore, a no-slip condition is assumed at the interface. 
The assumption of local thermodynamic equilibrium at the interface allows 
us the use of Raoult's law, in which the temperature dependence of the 
saturation pressure is determined via the Clausius-Clapeyron relation. 

The results present the temperature in the liquid and gas phase as well as 
the fluid motion in the liquid (caused by the evolution of the temperature via 
surface-tension and buoyancy effects) by means of streamlines as a function 
of time. The real total elapsed time is 10 seconds. Two videos are shown, 
presenting the same results in the following URLs: 



1. Video 1 - High resolution 



2. IVideo 2 - Low resolution! 

Note that in the videos the inner streamlines represent the highest velocity 
values. The red color represents the highest observed temperature (that 
of the ambient one), 298 K. The blue color represents the lowest observed 
temperature, around 285 K. 
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Figure 2: The temperature and liquid flow pattern at time t — Is for liquid 
thicknesses of 2 mm (left), 4 mm (middle) and 8 mm (right) and a gas flow 
of 100 ml/min 




Figure 3: The temperature and liquid flow pattern at time t — 10s for liquid 
thicknesses of 2 mm (left), 4 mm (middle) and 8 mm (right) and a gas flow 
of 100 ml/min 

2 Discussion 

From the results in the videos we can observe that first several small rolls 
are formed near the surface, caused by the surface-tension effect as Fig. [2] 
shows for the three liquid layer thicknesses at time t — Is. 
Due to buoyancy and as time proceeds, the rolls grow towards the bottom of 
the liquid layer. Then the rolls also grow in horizontal direction merging with 
each other until a steady configuration is obtained. For a higher liquid layer 
thickness, the merging occurs earlier and less rolls are left. Furthermore, 
the temperature gradients decrease as the liquid thickness increases, which 
is caused by the higher mixing efficiency when the liquid is less confined. 
Moreover, the rolls extend more horizontally under the cover towards the 
side walls as the liquid layer thickness increases. For smaller liquid layer 
thicknesses, the rolls reach the bottom where a constant temperature of 298 
K is maintained. Therefore the rolls stay concentrated close to the interface. 
As the liquid layer thickness increases, the rolls have more time to increase 
in size towards the side walls before they reach the bottom of the liquid 
layer. Fig. [3] shows this at the time t = 10s. 

This work yields valuable information about the supercritical instability 
behavior of an evaporating liquid and the qualitative influence of its con- 
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finement by means of fluid dynamics. 
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